UT-966 

Leptogenesis via LH U flat direction 
with a gauged U(1)b-l 



Masaaki Fujii, K. Hamaguchi 
Department of Physics, University of Tokyo, Tokyo 113-0033, Japan 



T. Yanagida 

Department of Physics, University of Tokyo, Tokyo 113-0033, Japan 

and 

Research Center for the Early Universe, University of Tokyo, Tokyo, 113-0033, Japan 



Abstract 



We study the supersymmetric leptogenesis via LH U flat direction with a 
gauged U(1)b~l symmetry. We find that the resultant baryon asymmetry is 
enhanced compared with the case without a gauged U(1)b~l symmetry. The 
baryon asymmetry is proportional to the reheating temperature of inflation, 
but it is independent of the gravitino mass. If high reheating temperatures 
of inflation Tr ~ 10 10 GeV are available, the mass of the lightest neutrino, 
m u i ~ 10 -4 eV, is small enough to explain the baryon asymmetry in the 
present universe. Furthermore, the gravitino mass independence of the pro- 
duced baryon asymmetry allows us to explain the present baryon asymmetry 
even in low energy SUSY breaking scenarios such as gauge-mediation models. 
Our leptogenesis scenario is also very special in the sense that it is completely 
free from the serious Q-ball formation problem. 
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I. INTRODUCTION 



Our universe is composed of baryon, but not anti-baryon. Accounting for this baryon 
asymmetry in the present universe is one of the most fundamental and challenging problems 
in particle physics as well as in cosmology. Various mechanisms have been proposed so far 
to solve this problem. In particular, there has been a growing interest in leptogenesis [0, 
since there are now strong evidences of neutrino oscillations, equivalently the existence of tiny 
but nonzero neutrino masses. It is, then, very natural to consider lepton-number violating 
interactions in nature, which is a crucial ingredient for leptogenesis to work. Especially, 
once we accept the supersymmetry (SUSY) and the existence of tiny Majorana neutrino 
masses, there automatically exist in the minimal supersymmetric standard model (MSSM) 
most of the necessary ingredients for the SUSY leptogenesis via LH U flat direction, which 
was originally suggested by Murayama and one of the authors (T.Y.) based on the idea 
of the Affleck-Dine (AD) mechanism [KJ. 

Recently, we performed a detailed analysis of this MY leptogenesis via LH U flat direction 
in the context of gravity-mediated SUSY breaking scenarios including all of the relevant ther- 
mal effects 0,|7|. We found a very interesting aspect of this leptogenesis, that is, "reheating- 
temperature independence of cosmological baryon asymmetry" J7|. The resultant baryon 
asymmetry is almost determined by the mass of the lightest neutrino (and the gravitino 
mass), and hence we can determine the lightest neutrino mass from the baryon asymmetry 
in the present universe, quite independently of the reheating temperature of inflation. This 
results in a crucial prediction on the rate of the neutrinoless double-beta (Oz//3/3) decay [0, 
which may be testable in future Qvf3[3 decay experiments such as GENIUS [H]- 

In this article, we study the effects of gauging the U{1)b-l symmetry on the leptoge- 
nesis via the LH U flat direction. (B and L are baryon and lepton number, respectively.) 
As well known, the U(1)b-l symmetry is the unique global symmetry which can be gauged 
consistently with the MSSM. If we accept a gauged U{1)b-l symmetry, the existence of 
three families of the right-handed neutrinos are automatically required from the anomaly 
cancellation conditions. Furthermore, if we assume that this U(1)b-l symmetry is sponta- 
neously broken at high energy scale, these right-handed neutrinos acquire large Majorana 
masses of the order of the B — L breaking scale, which naturally explains the tiny neutrino 
masses suggested from the recent neutrino oscillation experiments via the so-called "seesaw 



1 See also later studies 
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mechanism" ||. 

Below the B — L breaking scale, the scalar potential along the LH U direction is almost 
the same as in the case without a gauged U(1)b-l symmetry. However, above the B — L 
breaking scale, the scalar potential is lifted by the effect of the U(1)b-l D-term in a certain 
region of parameter space, which provides us a new mechanism to stop the <ft field. (</> 
is the superfield parameterizing the LH U "effectively flat" direction.) We show that if the <p 
field is stopped by the D-term potential, the baryon asymmetry is enhanced compared with 
the "global" case without a gauged U(1)b-l symmetry. This is because the available A-term, 
which provides a phase rotational motion of the field generating the lepton asymmetry, is 
much bigger than that in the global case. On the other hand, if the field is stopped by 
the F-term potential coming from the operator responsible for generating the small neutrino 
mass, the obtained baryon asymmetry becomes the same as that in the global case |||7]]. (In 
the following discussion, we will call the former case "D-term stopping case" , while the latter 
one "F-term stopping case.") 

Surprisingly enough, as we will see in the following discussion, if the amplitude of the (j) 
field is fixed by the D-term potential during inflation, the baryon asymmetry linearly depends 
on the reheating temperatures of inflation, but it is completely independent of the gravitino 
mass. This is totally an opposite situation to the global case, in which the final baryon 
asymmetry is almost independent of the reheating temperatures of inflation and linearly 
depends on the gravitino mass. 

We find in the D-term stopping case that if high reheating temperatures, Tr ~ 10 10 GeV, 
are available, the lightest neutrino of mass, m v i ~ 10~ 4 eV, can account for the present 
baryon asymmetry. Furthermore, the "gravitino-mass independence of the cosmological 
baryon asymmetry" provides us a great advantage in gauge-mediated SUSY breaking scenar- 
ios [JTTJ] . Gauge-mediated SUSY breaking is widely considered as one of the most interesting 
mediation mechanisms of SUSY breaking, since it can naturally explain the suppression of 
the flavor changing neutral currents (FCNC). However, in those scenarios, the reheating 
temperature of inflation is strongly restricted from above in order to avoid the overpro- 
duction of gravitinos [Q. This makes most of the baryo/leptogenesis scenarios impossible, 
especially when the gravitino is lighter than O(MeV), since in this case the reheating tern- 
perature T R is restricted as T R < 10 5 GeV.Q It was considered that the AD baryogenesis H 



2 The electroweak baryogenesis seems very difficult to take place in gauge-mediated SUSY breaking 
models because it is hard to obtain the required light stop, < rat- 
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(which is the baryogenesis based on the AD mechanism) is the unique candidate to produce 
an enough baryon asymmetry in gauge-mediated SUSY breaking scenarios. However, from 
a recent detailed analysis, it becomes clear that the AD baryogenesis is not successful in 



gauge- mediated SUSY breaking scenarios because of a serious Q-ball formation fll3[ . In the 
present work, we will show that the MY leptogenesis can naturally explain, in the D-term 
stopping case, the present baryon asymmetry even in gauge-mediated SUSY breaking mod- 
els with the gravitino mass smaller than (9(100 keV), if the mass of the lightest neutrino is 
m u i ^ 10~ 10 eV. We should stress that such a light gravitino may be directly detected in 
future collider experiments 0]. We claim that our scenario is the first minimal model which 
explains the present baryon asymmetry in the gauge mediation model with a small gravitino 
mass m 3 / 2 = C(100 keV). 

II. THE MODEL 

First, let us discuss the potential for the relevant fields. To demonstrate our point, we 
use a simple superpotential as follows: 

W = fiH u H d + hNLH u + -XSNN + 7]X(SS-v 2 ) . (1) 

2 

Here, h, A, 7] are coupling constants and we assume rj = 0(1). X, S and S are singlets under 
the MSSM gauge group and S, S carry —2, +2 of B — L charges, respectively. N, L, H u 
and Hd are a right-handed Majorana neutrino, an SU (2) L-dovh\et lepton and Higgs fields 
which couple to up and down type quarks, respectively, v is the breaking scale of the U(l) b-l 
symmetry.^ As we will see later, the "effectively flat" direction relevant for the present baryon 
asymmetry is the flattest flat direction, and hence we will consider only one family which 
corresponds to the lightest left-handed neutrino. We adopt the following linear combination 
of L and H u as the flat direction field <fi : 

/n \ (2) 

The Kahler potential must have non-minimal couplings of the <fi field to the inflaton, 
otherwise the <fi field gets a large positive mass term of the order of the Hubble parameter 





3 Here and hereafter, we take the couplings h, A, rj and /i, v to be real. This can be done by field 
redefinitions in Eq.(ffl). 
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(which we will call a Hubble mass term) and it is driven exponentially towards the origin 
during inflation, and hence the leptogenesis via LH U flat direction cannot take place ||. 
Therefore, we assume that there are general non-minimal couplings of the <fi field to the 
inflaton in the Kahler potential, 

SK 3 (jfclfo + h.c.) + + • • • • (3) 

Here \c^\, = 0(1) are coupling constants, and M* ~ 2.4 x 10 18 GeV is the reduced Planck 
scale. / is the inflaton superfield. From these non-minimal couplings to the inflaton, the </> 
field gets the following SUSY breaking effects during inflation and at the inflaton-oscillation 
dominated epoch: 

5V 3 -V3#MW* + h.c.) + 3(1 - b'^H 2 ^ 2 + . . . , (4) 

where H is the Hubble parameter, b'^ = b^ — |c^| 2 , W$ = dW/d(p, and we use the fact that 
|F/| 2 ~ 3H 2 M 2 in these regimes.n In the following discussion, we assume 3(1 — b'A ~ — 1 
for simplicity. We also assume that there are non-minimal couplings of other fields to the 
inflaton in the Kahler potential. Then, the full scalar potential relevant to the field is given 

byQ 



V = -g 2 (-2|S| 2 + 2|£| 2 + |Af -|L| 2 ' 2 

|2 , I 7 T ZT , \ O T\T\2 



+\ptH u \ z + \hLH u + XSN\ + \fiH d + hNL\ + \hNH u \ 2 
+ ±\N 2 + V XS + \rj(SS -v 2 )\ 2 + \nXS\ 2 

+3H 2 (^(l-b'y)^ -VSH ^CyYWy + h.C^j+VsB 



(5) 



where b' Y = by — \cy\ 2 , cy are non-minimal couplings of Y(= L, H u , Hj, X, S, S, N) to the 
inflaton in the Kahler potential as coupling constants in Eq. (H), and g = 0(1) is the gauge 
coupling constant of the U(1)b-l- Vsb represents soft SUSY breaking terms in the present 
true vacuum. Here, all fields denote the scalar components of the corresponding superfields. 



We redefined the phase of c^. 



5 Here, we omit the potential coming from thermal effects, which will be discussed in Section IV 
and[V|. 
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The scalar potential in Eq. (|5|) is so complicated that it seems very difficult to solve the 
dynamics of the relevant fields. Fortunately, however, we only need to know the shape in 
the neighborhood of the bottom of this potential, if the curvature around this bottom is as 
large as the Hubble parameter during inflation. This is because all the scalar fields, which 
have masses as large as the Hubble parameter, settle down at the bottom of the potential 
during inflation and trace this potential minimum throughout the history of the universe. 
This allows us to eliminate many terms in the potential by using the field relations required 
for minimizing the potential. Therefore, let us first find out the conditions to minimize the 
potential in Eq. (||). Though this is also a very hard task because of the complexity of the 
potential, one can find out the minimum of the potential at least when a condition |0|, H^zv 
is satisfied. In fact, we find the approximate minimum of F-terms for H ^ v, [] 

S^j, X^O, N~-^ } H d ^0. (6) 

The curvatures around the first three minimums are of the order of the B — L breaking scale 
v, and that around the last one is of the order of the Hubble parameter. As explained above, 
we simplify the potential by using the relations in Eq (||). Then, the potential in Eq. ([5]) is 
reduced to 

?; 4 h 
-2|S| 2 + 2 ^ + 





S\ 2 AX 2 \S\ 2 

-H 2 ^ 2 + 3(1 - b' s )H 2 \S\ 2 + 3(1 - b' B )H 2 ^- + 3(1 - b' N )H 2 h * 



\S\ 2 y NJ A\ 2 \S\ 
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+ ^(c;y^ + h.c.)+^ 2 |0| 2 + ^. (7) 

Here, c'^ = — 1/4 cs- One may wonder why the "flat" direction field can develop a large 
expectation value in spite of the presence of the U{1)b-l D-term. This is because the S field 
shifts and absorbs the D-term potential. Here, we have 



6 Here, we assume a positive Hubble-order mass term for the Hd field. This is necessary in order 
to avoid contamination of to the LH U flat direction. We have checked both analytically and 
numerically that the contamination of the field becomes nonnegligible only after the H Jo and 
hence the following discussion is not affected much by the contamination. 
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< ! li + ^W- (8) 



Then, we get the following potential: 



V ^ V SB + ^ 2 |0| 2 + ^ ( </ * £ + h.c.l + * £L + O 



-tf 2 |0| 2 + 3(1 - b' s )H*\S\ 2 + 3(1 - ^)# 2 ^ + 3(1 - b' N )H^J-^ 2 . (9) 

The above procedure can be justified as long as the B — L breaking scale v is larger than the 
Hubble parameter during inflation, because the curvature around the minimum in Eq. (|8|) 
coming from the D-term potential is also of the order of the scale v. We assume that this 
is the case in the following discussion, which is, however, very plausible since the B — L 
breaking scale v and the Hubble parameter of inflation Hj are most likely v ^ 10 13 GeV and 
Hj ^ 10 13 GeV. As we will see later, the third term, 

plays a crucial role in the D-term stopping case. We will call this term "Hubble A-term" in 
the following sections. 



III. EVOLUTION OF THE FIELD DURING INFLATION 

In this section, we discuss the evolution of the field during the inflation. As we will see 
later, in the D-term stopping case, the field settles down at the B — L breaking scale v 
during the inflation. This allows us to use the Hubble A-term as a torque to rotate the field 
after the inflation ends. It will become clear in Section ^ that the torque coming from this 
Hubble A-term is the crucial ingredient to enhance the baryon asymmetry compared with 
the global case as well as the F-term stopping case. 

First, let us investigate the potential for the field during the inflation in the range 
|0| ^ v. In this range, the S field in Eq. (§) is expanded into the following form: 

|S| 2 ^ 2 -i|0| 2 + o(^, for \<f>\<v. (11) 

Then, by substituting Eq. (|TTf) into Eq. (||), we get the effective potential for the field in 
the range |0| in the following form: 
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1 , IJL , a , VSH f , ^ x |0| 



^-^ + ^ 2 i0i 2 + ^^(c;0 4 +h.c.) + 



(i 



2 a- in 2 M V ^ ' V 4M 2 

-^ 2 i0i 2 +^(^-^)^vr+-- - , (12) 

where the ellipsis denotes higher order terms in |0| 2 /t> 2 . Here, M = Xv/h 2 corresponds to 
the same symbol used in Ref. [|],||[7| and the mass of the lightest neutrino is written as 
m ul = (H u ) 2 /M. 

From Eq. (|12|), we see that the effective potential is almost the same as in the global case 
with some changes of coefficients in the potential. However, there is an important difference 
in the Hubble-induced SUSY breaking mass terms. We should not neglect the contributions 
coming from the S field. In fact, we need the following condition in order that the field 
develops a large expectation value during the inflation: 

b' S -b's<l- (13) 

If this is the case, the field has a negative Hubble mass term at least in the range |0| ^ v. 
Otherwise, the field is driven toward the origin during the inflation because of the positive 
Hubble mass term, and the MY leptogenesis does not take place. 

In the following discussion, we assume that Eq. flTB| ) is satisfied and discuss the scale 
where the flat direction is lifted. If the balance point between the negative Hubble mass term 
and the F-term potential |0| 6 /(4M 2 ) in Eq. ( |T2D is below the B — L breaking scale v, the 
field is stopped at this balance point, 



'MHi < v , (14) 
which corresponds to the F-term stopping case. On the other hand, if the condition 



MHj>v, (15) 

is satisfied, the field can develop its expectation value as large as v. 

Now, let us explain what happens when the expectation value of the field grows as large 
as the B — L breaking scale v. At this scale {i.e. |0| ~ v ), the expansion of the S field given 
in Eq. QTTj ) becomes invalid and above this scale we must use another expansion of the S 
field as follows: 



|S| 2 -4^ + |_|0| 2 + 




for |0| >v. (16) 
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One might wonder whether the above expansion is reliable, since Eq. (|j) is based on Eq. (^) 
that may not be applicable for \cf)\ ^> v. Here, we first derive the conditions to fix the <fi 
field at the scale v during the inflation, assuming this expansion is effectively applicable at 
least for |0| ~ v. We will justify later the validity of the obtained conditions by numerical 
calculations. 

By substituting the expansion Eq. ([TJ]) into Eq. @, we obtain the following Hubble- 
induced mass term:[] 

V^^-l + 3 -(l-b^H 2 \^ + o(^H^ for \<j>\>v. (17) 

Then, the <fi field gets a positive Hubble mass term for \<fi\ ^ v if the following condition is 
satisfied: 



Therefore, if this is the case, the <\> field cannot develop its expectation value above the scale 
v and hence it is fixed at the B — L breaking scale v during the inflation. 

To summarize, the <fi field, and hence S and S fields also, are stopped at the B — L 
breaking scale v during the inflation, if the following conditions are satisfied,^ 

Ks-b's<l JmHi>v, &'s<-~" ( 19 ) 

These are the conditions for the D-term stopping case. 

To check the conditions in Eq. (|l9|), (especially that of the last condition b' s ^ —1/3,) we 
have numerically solved the coupled equations of motions for the relevant fields using the 
full scalar potential in Eq. (^). We show the result in Fig. [I], where the amplitude of the (j) 
field at the end of the inflation is plotted in b' s -b'g plane. It is found that the conditions in 
Eq. ( |T9D well explain the result of this numerical calculation. In fact, the amplitude of the 4> 
field at the end of the inflation lies in the range l^S |0|/f^$3 in most of the parameter space 
where the conditions in Eq. flTPp are satisfied. 



7 Here, we assume h 2 /X 2 < 1. If this is not the case, we must include the Hubble mass term 
coming from the coupling of the right-handed Majorana neutrino to the infiaton, which is the last 
term in Eq. (H|). 

8 The conditions for the region h 2 /X 2 > 1 are obtained by repeating the same procedure. 
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As noted in the first paragraph of this section, we can use the Hubble A-term potential 
[Eq. ([lOD l as a torque to rotate the (f> field after the inflation ends in the D-term stopping case. 
In the remaining part of this section, we discuss this point in detail and compare it with the 
F-term stopping case. 

First, if the conditions in Eq. ([19]) are satisfied {i.e., the D-term stopping case), the 
curvature of the potential for the field along the phase direction (which is denoted by the 
symbol mp tase ) is smaller than the Hubble parameter during the inflation. This can be seen 
from the following relation: 



-4- - >' = ( JL) 2 < H] . (20) 

Then, there is no reason to expect that the <ft field sits down at the bottom of the valley of 
the A-term potential in Eq. (10) during the inflation. Thus, unless there is an accidental 
fine tuning on the initial phase of the <fi field, it is generally displaced from the bottom of 
the valley of this A-term potential when the inflation ends. Therefore, this A-term potential 
kicks the (f) field along the phase direction when the Hubble parameter becomes comparable 
with the curvature along the phase direction mp hase . This is the reason why we can use the 
Hubble A-term potential as a torque to rotate the <fi field and enhance the baryon asymmetry. 

Next, we turn to the F-term stopping case. In this case, the amplitude of the field 
is fixed at ^MHi during the inflation. [See Eq. (0).] An important point is that the 
curvature around the valley of the Hubble A-term potential, mL se , is of the order of the 
Hubble parameter with the value of the <p field, \4>\ — \fMHj\ 

Mphase ^^F\4> 2 \^H} . (21) 



Therefore, the phase of the <fi field settles down at the bottom of the valley of this A-term 
potential during the inflation, and hence the Hubble A-term cannot supply a torque to rotate 
the (j) field. In this case, the relevant torque for the field only comes from the ordinary 
A-term potential proportional to the gravitino mass|] In this case, the MY leptogenesis with 
a gauged U(1)b-l symmetry results in the same conclusions as in the global case H0], in 
which the resultant baryon asymmetry is proportional to the gravitino mass. We stress that 



9 The initial phase of the <p field is generally displaced from the bottom of the valley of this ordinary 
A-term potential unless the valley of this A-term potential accidentally coincides with that of the 
Hubble A-term potential. 
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this leads to a substantial suppression of the baryon asymmetry in gauge-mediation models 
with a small gravitino mass.0 

Finally, we comment on the case where only the condition b'g < —1/3 in Eq. ( fT9f) is 
not satisfied. In this case, the <fi field has the negative Hubble-mass term even above the 
B — L breaking scale. We need the effective potential for the </> field above the scale v in 
order to determine the scale at which the (ft field stops during the inflation. However, there 
are so many coupling constants and terms in Eq. (|5]) which contribute to the potential with 
comparable importance and hence it is very hard to get a simple effective potential. In spite 
of this complexity of the potential, it is expected that the baryon asymmetry results in the 
range between the following two cases. If the curvature around the valley of the Hubble 
A-term potential is as large as the Hubble parameter during the inflation, mp kse ~ Hj, the 
baryon asymmetry approaches what is obtained in the global case [§,0. On the other hand, 
if mp hase < Hj, the baryon asymmetry is enhanced and probably the same as that in the 
D-term stopping case with which we are mainly concerned in this paper. In any case, we need 
detailed numerical calculations using the full scalar potential in Eq. (||) in order to obtain 
the precise amount of the baryon asymmetry in this case, which is beyond the scope of this 
work. 

IV. BARYON ASYMMETRY 

In this section, we calculate the baryon asymmetry. First, let us discuss the D-term 
stopping case. In this case, as explained in the previous section, the <p field is fixed at the 
B — L breaking scale v during the inflation, where the curvature along the phase direction 
coming from the Hubble A-term is smaller than the Hubble parameter Hi. Thus, the phase 
of the field is generally displaced from the bottom of the valley of the Hubble A-term 
potential. 

After the inflation ends, the universe becomes dominated by the oscillating inflaton, and 
the scale factor of the expanding universe R increases as R oc H~ 2 ^ 3 like in the matter- 
dominated universe ]18[. In the following discussion, we assume that the production of the 
lepton number takes place in this inflaton-oscillation dominated epoch, before the reheating 
process of the inflation completes. (We will justify this assumption in the next section.) 



See the discussion in the next section. 
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At first, the field value of the remains almost constant for a while. This is because 
the amplitude of the held which gives the potential minimum is mainly determined by the 
second line in Eq. (H) and it is given by |0| ~ v which is, of course, independent of the Hubble 
parameter H. 

Then, the Hubble A-term potential in Eq. ( |10|) kicks the field along the phase direction 
when the Hubble parameter becomes comparable with the curvature. The Hubble parameter 
of this time is given by 

Bv^Jj # Hl h * mj hase ^fV) . (22) 

At this time (H ~ H p h), the phase of the field begins to oscillate around the bottom of 
the valley of the Hubble A-term potential. On the other hand, the amplitude of the field 
slowly decreases after this time, according to 

\<j)\~VMH, (23) 

which is the balance point between the negative Hubble mass term and the operator 
|0I7(4M 2 ). 

Since the field starts its oscillation along the phase direction at H = H p h, it has already 
had a large acceleration along the phase direction before it starts coherent oscillation along 
the radius direction around the origin when H = H osc . This H osc is determined by thermal 
effects, soft SUSY breaking effects and /i-term in the same way as in Ref. ]/]]. For example, if 
the field starts its oscillation due to the soft SUSY breaking mass in the true vacuum m^, 
it is given by H osc ~ m^. The evolution of the amplitude of the field and its phase obtained 
by numerical calculations are shown in Fig. |2] and || In these figures, t is the cosmic time in 
the matter dominated universe t = 2/(3H). We see from these figures that the evolution of 
the field is well explained by above arguments. 

We are now at the point to estimate the baryon asymmetry. The lepton number density 
is related to the field as follows: 

n L = ^(0*0-0*0) , (24) 

where the overdot denotes a derivative with time. The evolution of the field is described 
by the following equation of motion: 

+ 3^0+|^ = O. (25) 
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Then, we obtain, from Eq. (]24]) and (^), the equation of motion for the lepton number 
density as0 

n L + 3Hn L = ^Vs^lm (c^ 4 ) . (26) 

Here, we have used \S\ ~ v for H^z H ph . This can be easily integrated and we obtain the 
lepton number at time t as 

(*) = /* dtR 3 2 V / 3^V| sin (arg( C ;) + 4arg(0)) . (27) 

In the regime H osc ^ H ^ H phl the amplitude of the field decreases as \<fi\ ~ \/MH and the 
sign of sin in Eq. ( |2T| ) changes with time scale. (See Fig. |2] and |3|.) Therefore, from 
Eq. (^), one sees that the total lepton number oscillates with 1/H time scale with almost 
constant amplitude, since i? 3 if |0 4 | oc H oc t~ l in this regime. As soon as the field starts 
its coherent oscillation around the origin at H = H osc , the total lepton number is fixed, since 
the amplitude of the field decreases as fast as \<f>\ oc H. The lepton number density at time 
t = 2/(3H) is written as 

n L (t) * ^=\c^\MH% ff , (28) 

where 5 e ff = sin (arg(c^) + 4arg(0)) is an 0(1) effective CP-phase. Note that the Eq. ( |28|) 
is correct as long as H ^ H ph , and hence the resultant lepton-to-entropy ratio is totally 
independent of the time when the field starts oscillations around the originQ 

V± „ ^Lu\ S „ (29) 
where Tr is the reheating temperature of the inflation. This lepton asymmetry is partially 



converted into the baryon asymmetry |l| due to the "sphaleron" effects 15 , since it is 



11 Here, we neglect the contribution coming from the ordinary A-term, since its contribution is 
negligible. 

12 We have also confirmed by numerical calculations that the final lepton asymmetry n^/s is inde- 
pendent of the oscillation time H osc . See also the result for the F-term stopping case in Eq. (|34]). 
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produced before the electroweak phase transition. The present baryon asymmetry is given 
by 



(30) 



n B _ 8 n L 
s ~ 23 s 

Thus, after all, the present baryon asymmetry is given by 

ub 8 T R M 
s ~ 23VS Ml |C ^ e// 

Here, we have used the relation M = (H u ) 2 /m v \. 

As stressed before, one sees that the resultant baryon asymmetry is independent of the 
gravitino mass 7713/2 and the starting time of the oscillation of the field, H osc . If the mass 
of the gravitino is large enough m 3 / 2 — 0(1 TeV), we can avoid the cosmological gravitino 
problem even if the reheating temperature is rather high Tr ~ 10 10 GeV |I7|]. In such a 



case, we see from the above equation that the mass of the lightest neutrino m u i — 10~ 4 eV 
is small enough to explain the present baryon asymmetry On the other hand, if the mass 
of the lightest neutrino is as small as m v \ ~ 10~ 10 eV, we can generate the required baryon 
asymmetry for low reheating temperatures as Tr ~ 10 4 GeV. In such a low reheating 
temperature, we are free from the overproduction of gravitinos even in gauge-mediated SUSY 
breaking scenarios with the small gravitino mass 777.3/2 = 0(100 keV) |T2"f . 



We show the evolution of the baryon asymmetry — = obtained by a numerical 

s 23 s 

calculation in Fig. [| and compare it with the estimated value in Eq. fl3"T|). From this figure 
we see that the asymmetry begins to oscillate with almost constant amplitude at H = H p h 
and it is fixed at the time when the field starts to oscillate around the origin. We confirm 
from this figure that the resultant asymmetry is well explained by the arguments described 
in this section. 

Next, we briefly discuss the F-term stopping case, in which both the evolution of the <fi 
field and the resultant baryon asymmetry are completely the same as those in the global 
case R. In this case, the (f) field starts to move as soon as the inflation ends, according to 



13 In the present analysis, we neglect the relative sign between the produced lepton and baryon 
asymmetries. 
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|0| ~ V MH . As mentioned in the previous section, the Hubble-induced A-term cannot play 
a role to kick the phase of <f), since the flat direction field <fi is trapped in the valley of this 
Hubble-induced A-term in this case. The only available A-term that can cause the motion 
of cf) along the phase direction is the ordinary A-term, which is proportional to the gravitino 
mass m 3 / 2 : 

VsB3 ^m ( a ^ 4 + h - c -) > ( 32 ) 

where a m is an 0(1) constant. Thus, the equation of motion for the lepton number density 
in Eq. fl26|) is replaced by 

h L + 3Hn L = 7^-I m ( a ™0 4 ) • (33) 

By repeating the same procedure as in the D-term stopping case, we obtain the following 
baryon asymmetry: 

2 T R M fm 3/2 



n B 
s 



F-term 69 M? V H 



\a m \8eff • (34) 



The factor (m 3 / 2 /H osc ) gives rise to a strong suppression of the resultant baryon asymmetry 
for high reheating temperatures, since H osc becomes much larger than the soft mass in 
that region |||7[].[3 We also stress here that this factor shows the strong suppression of the 
baryon asymmetry in the gauge-mediated SUSY breaking scenario, since m 3 / 2 <C < H osc . 



V. OTHER CONSTRAINTS 



In the previous section, we found that the resultant baryon asymmetry is enhanced in 
the D-term stopping case compared with the F-term stopping case and also the global case. 
However, in the D-term stopping case, there are some conditions we must check, other than 
those in Eq. (|H5|). These are the conditions for avoiding the early oscillation of the <fi field. 
Although the baryon asymmetry in the D-term stopping case is independent of the H osc as 
long as H osc < H p h, if the <fi field starts its oscillation around the origin before H = H p h, 



The "reheating temperature independence of the baryon asymmetry" in the global case comes 
from the fact that H osc in Eq.(|34|) is proportional to Tr or T^ 2 in a large part of the parameter 
space, which leads to ns/s oc or T^ 3 Q. 
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i.e. H osc > H p h, the resultant baryon asymmetry is strongly suppressed.^ The soft SUSY 
breaking mass term of the field and the //-term cause the early oscillation of the field 
when they satisfy^ 

v 2 

m h fj, > H ph ~ — . (35) 

Furthermore, thermal effects may cause the early oscillations. A field which couples to 
the field gets an effective mass fk\4>\, where fj. is a coupling constant to the field. If the 
cosmic temperature T is larger than this effective mass, thermal fluctuations of that field 
produce the thermal mass term for the field, Cfc/|T 2 |0| 2 . The list of c& for the LH U flat 
direction is given in Ref. ||. Then, if one of these thermal mass terms exceeds the Hubble 
parameter in the regime H > H p h, it causes the early oscillations of the field, and then 
suppresses the baryon asymmetry substantially. Namely, the following condition is needed 
to avoid the early oscillation of the field by the thermal mass terms: 

H 2 > E c kfl T2 for H>H ph . (36) 

f k v<T 



There is another thermal effect which was pointed out in Ref. ||19|| . Along the LH U 
flat direction, SU (3) gauge symmetry remains unbroken, and hence, gluons and gluinos are 
massless. Furthermore, the down type (s) quarks also remain massless since they have no 
coupling to the field. These light fields produce the free energy which depends on the 
S77(3) gauge coupling constant. We obtain the effective potential V cx #§T 4 at two loop 
level |[20|| , where gs is the gauge coupling constant of the 577(3). At first sight, there seems 
no dependence on the field in this free energy. However, the up type (s) quarks get large 
masses from the couplings to the field y u , and if y u \4>\ > T, they decouple and change the 
trajectory of the running coupling constant of the 577(3). This effect produces the effective 
potential for the field,Q 



15 In this case, the resultant baryon asymmetry is inversely proportional to H 2 SC . 

16 Note that, in gauge-mediated SUSY breaking scenarios, the condition for is absent if the 
messenger scale is lower than the breaking scale of the U(1)b-l- The existence of the /i-term 
depends on the energy scale of the dynamics which produces the /i-term. 

17 For non-abelian gauge symmetries, decouplings of matter particles give always positive contri- 
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5V ~ a,a|T 4 log ^ , (37) 



where a g is a constant a bit larger than unity and as = 5 , |/4vr J7|,P^[. Then, the following 
condition is required to avoid the early oscillation caused by this potential: 



2 rpA 

H 2 > Qga f for H> H ph . (31 



From Eq.fl36|) and (p8|), we obtain the following condition to avoid the early oscillation by 
thermal effects: 



Tr < min 



minfc < max 



(39) 



where, we have used the fact that the cosmic temperature behaves as T = (ifT^M*) 1 / 4 before 
the reheating process of the inflation ends [TJ[]. Taking all of these effects into account, 
we get the allowed regions which are free from the early oscillations as in Fig. |5|. The 
regions below four solid lines are free from the early oscillations. These four lines correspond 
to the breaking scale of the U(1) B -l, v = 10 16 , 10 15 , 10 14 , 10 13 GeV from left to right, 
respectively. The vertical parts of these lines come from the condition in Eq. fl35|). Here, 
we assume the existence of or /x-term when the amplitude of the <fi field is of the order 
of the B — L breaking scale. The shaded region denotes the present baryon asymmetry, 
Ub/ 's ~ (0.4 — 1) x 10~ 10 . From the Fig. |5|, we see that the early oscillations can be easily 
avoided, especially when the B — L breaking scale satisfies v ^ 10 14 GeV.0 

Finally, we here check the assumption that the production of the lepton asymmetry takes 
place during the inflaton-oscillation dominated era, i.e., before the reheating process of the 
inflation completes. Then, we have the following constraint; 

H >r,~f^Y /2 ^ (40) 



butions, and the net contribution from decouplings of gauge bosons and gauginos is always negative. 
The LH U flat direction has only positive contributions for SU(3) gauge symmetry. On the other 
hand, other flat directions such as udd, LLe, . . . receive large negative contributions and probably 
have negative thermal-log term. 

18 In the regions above the four solid lines, the estimation of H osc is different from that in Ref. Q. 
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where Tj is the decay rate of the inflaton / and g* denotes the number of relativistic degrees 
of freedom, which is g*{T) ~ 200 in the MSSM for T ^> 1 TeV. As long as there exists a 
soft mass or the //-term at high energy scale, H osc is at least larger than that. Thus, the 
above constraint is indeed satisfied for T R ^ 10 10 GeV x (m^or/z/l TeV) 1//2 . Furthermore, 
the Hubble parameter of the oscillation time H osc becomes much larger for higher reheating 
temperatures, because of the thermal effects. For example, the H osc always satisfies the 
following relation 0; 

/ M \ 1/2 

Hosc > a s T R [a^j . (41) 
Thus, the constraint given in Eq. (fl0|) is satisfied as long as 

T " ~ 1015 GeVx (upFw) • < 42 > 

which is the case in all the relevant parameter space in the present analysis. (See Fig. |5|.) 
Therefore, the assumption Tj < H osc is justified. 



VI. Q-BALL PROBLEM 

One of the physical reasons to make the leptogenesis via LH U flat direction so special 
among other baryo/leptogenesis scenarios using the AD mechanism is that it is free from 
the Q(L)-ball |21[] problem. The coherent oscillation of the <fi field is unstable with spatial 
perturbations if the potential of the field is flatter than the quadratic potential.^ If this 
is the case, the coherent oscillation of the field fragments into non-topological solitons, Q- 
balls ||22|| . By the recent detailed analysis by lattice simulations, it becomes clear that almost 
all of the charges carried by the field are absorbed into these Q-balls p3| , and the present 



baryon asymmetry must be provided by the decay of the Q-balls, not by the direct decay 
of the (f) field. The formation of Q-balls is a generic feature of the AD baryo/leptogenesis 
regardless of the mediation mechanism of SUSY breaking. 

In gravity-mediated SUSY breaking scenarios, the potential of the field is slightly flatter 
than the quadratic potential due to the running of the soft mass of the field coming from 



19 In this section, we will use the same symbol <fi to denote the field which parameterizes a general 
flat direction. 
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gaugino loops |24|,^]. In this case, the Q-ball is generally unstable. A difficulty arises from its 
long life time. If the decay temperature of the Q-balls is well below the freeze out temperature 
of the lightest supersymmetric particle (LSP), the amount of the LSP cold dark matter can 
be written as 

°«- 3 (43) 

where N x is the number of LSP's produced per baryon number, which is at least 3, and 
fs — 1 is the fraction of baryon number stored in the form of Q-balls, and m n and m x are 
the nucleon mass and the neutralino (%) LSP mass, respectively. Qx denotes the ratio of 
the energy density of X to the critical density of the present universe. From this relation, 
it is clear that the late time Q-ball decay leads to the over production of the LSP cold dark 
matter. In fact, this late time decay of Q-balls invalidates most of the AD baryo/leptogenesis 
scenarios in gravity-mediated SUSY breaking models.0 

On the other hand, in gauge-mediated SUSY breaking scenarios, the soft SUSY breaking 
mass of the <fi field vanishes above the messenger scale, and the (f> field has only logarithmic 
potential. This fact leads to the formation of Q-balls. The Q-balls in gauge-mediated SUSY 
breaking scenarios have very different characters compared with those in gravity-mediated 
SUSY breaking models. The effective mass of the Q-ball per baryon number is proportional 
to Q -1 / 4 ||26||, where Q is the charge of the Q-ball. Therefore, if the charge of the produced 



Q-ball is large, the Q-balls become completely stable and they remain as a cold dark matter 
in the universe Once we fix the flat direction for the AD baryo/leptogenesis, we can 

estimate the size of the produced Q-ball and its number density by using the result of the 
recent numerical calculations in Ref. We find that the produced Q-balls overdose the 



universe in almost all of the regions of parameter space. Though there remains some tiny 
regions in parameter space in which the produced Q-balls do not overdose the universe, 
most part of such regions are already experimentally excluded [[OJ. If we use the leptonic 
flat directions, such as LLe, the produced Q(L)-balls can decay into neutrinos, and hence 
they do not overdose the universe. Unfortunately, however, we can only use the leptonic 
charges which evaporates before the electroweak phase transition, because we must convert 



20 The present status of AD baryo/leptogenesis in gravity-mediated SUSY breaking models is sum- 
marized in Ref. There, we also proposed an interesting AD baryo/leptogenesis model with 
a gauged U(1)b-l symmetry, which is an almost unique solution to solve the Q-ball problem in 
gravity-mediation models except the present leptogenesis via LH U flat direction. 
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the lepton asymmetry to the baryon asymmetry by "sphaleron effects" ||15|| . As a result, we 
find that the present baryon asymmetry can be explained only in extremely small regions of 
the parameter space .0 

How about the LH U flat direction? The soft SUSY breaking mass of the <fi field has no 
contribution from gluino loops, but, on the other hand, it has a big opposite contribution from 
large top Yukawa coupling, which makes the potential for the field steeper than quadratic 
potential. Therefore, the gravity-mediation type Q-balls are not formed in the LH U flat 
direction |[25|| . Furthermore, the LH U flat direction is the unique flat direction available for 
the AD baryo/leptogenesis which has a SUSY mass term, /i-term. Thus, once we assume 
the existence of the /i-term even in the high energy scale relevant for the MY leptogenesis, 
there always exists the quadratic potential for the field and the formation of Q-balls is 
avoided even in the gauge-mediated SUSY breaking scenarios. (This is possible if the /i-term 
is generated by an expectation value of a nonrenormalizable operator.) 

Therefore, the MY leptogenesis is free from the Q-ball problem not only in gravity- 
mediated SUSY breaking scenarios but also in gauge-mediated SUSY breaking scenarios]^] 



VII. DISCUSSION AND CONCLUSIONS 



In this paper, we have investigated the supersymmetric leptogenesis via LH U flat direc- 
tion originally suggested by Murayama and Yanagida (MY) 0] in the presence of a gauged 
U(1)b-l symmetry. We have shown that if the flat direction field <fi is stopped by the F-term 
potential, the situation is almost the same as that in the global case where the U{1)b-l is 
not gauged 0. 

Interestingly, however, we have found that the resultant baryon asymmetry in the D-term 
stopping case is enhanced compared with the global case. In this case, the field is stopped 
at the B — L breaking scale, where the curvature along the phase direction coming from the 
Hubble-induced A-term is smaller than the Hubble parameter. This allows us to use a much 



The authors thank S. Kasuya and M. Kawasaki for useful discussion. 

22 If the 4> field starts to oscillate around the origin due to the thermal-log term 5V = 
a 9 a 2 T 4 log ^|0| 2 /T 2 ^ , gauge-mediation type Q-balls are formed. Fortunately, however, in such a 
case, the Q-balls produced by the thermal-log term are small enough to evaporate well above the 
electroweak phase transition. Therefore, the thermal-log term does not harm the leptogenesis via 
LH U flat direction. 
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bigger torque to rotate the field, since the field is initially displaced from the bottom 
of the valley of the Hubble A-term potential unless we tune the initial phase of the (f) field. 
We have shown that the baryon asymmetry is proportional to the reheating temperature 
of the inflation, but it is completely independent of the size of the gravitino mass. As a 
result, if high reheating temperatures Tr ~ 10 10 GeV are available, the mass of the lightest 
neutrino, m v \ ~ 10 -4 GeV, is sufficiently small to explain the present baryon asymmetry. 
Furthermore, the gravitino mass independence of the resultant baryon asymmetry gives us 
a great advantage in gauge-mediated SUSY breaking scenarios. If the mass of the lightest 
neutrino is as small as m u i — 10~ 10 eV, we can explain the present baryon asymmetry 
even in the case where the gravitino mass is O(lOOkeV), which may be directly tested in 
future experiments ||14|1 . We stress that our scenario is the first minimal framework which 
can explain the present baryon asymmetry in gauge-mediated SUSY breaking scenarios with 
such a small gravitino mass. 

In addition, the MY leptogenesis via LH U flat direction is very special in the sense that it 
is completely free from the Q-ball problem. In gauge-mediation scenarios, this requires the 
existence of the /i-term (SUSY-invariant mass for Higgs multiplets) in the high energy scale 
relevant for the leptogenesis, which may imply the high energy origin of the /i-term in gauge- 
mediated SUSY breaking scenarios. The generation of the correct size of the /i-term has been 
one of the most difficult problems in gauge-mediated SUSY breaking models. The presence 
of the baryon asymmetry in our universe can be regarded as an interesting implication for 
the origin of the /i-term. 

Finally, we briefly comment on the prediction on the rate of the 0i>(3(3 decay. As in 
the global case, the present scenario also requires a hierarchical neutrino mass spectrum to 
generate the present baryon asymmetry even in the D-term stopping case. This results in 
almost the same prediction on the rate of the Qv(3(3 decay as that in the global case j?J. This 
consistency check of the present leptogenesis scenario will be given in future Qvf3f3 decay 
experiments ||. 
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FIGURES 
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FIG. 1. The amplitudes of the <ft field fixed during inflation which are determined by numer- 
ical calculations. In this calculation, we have used the full scalar potential in Eq. @ to follow 
the evolution of the relevant fields. Here, we have assumed that Hj/v = 0.1, 3(1 — b'^) = — 1, 
h = 10 -4 , g,X,rj = 1, and have randomly generated other coupling constants by, \cy\ in the range 
—2.5 < by < 2.5 and < \cy\ < 2.5, respectively. Various symbols denote the \4>\ to v ratio at the 
end of inflation. We see that the <p field is, in fact, stopped at the B — L breaking scale v if the 
conditions in Eq. (|l9|) are satisfied. 
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FIG. 2. The evolution of the <j) field estimated by a numerical calculation. We assume that 
M = 3 x 10 23 GeV (i.e., m vX = 1(T 10 eV), u = 10 14 GeV, Hj = 10 12 GeV, = 10 3 GeV. We see 
that the (f> field starts to move at H = H v h ~ 3 x 10 4 GeV according to the equation, \<$>\ ~ V MH . 
Here, t = 2/{3H) is the cosmic time. The (ft field starts to oscillate around the origin at H osc ~ 
since we neglect the thermal effects in this calculation. 
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l0g10 (ceN 

FIG. 3. The evolution of the phase of the cj> field estimated by a numerical calculation. The 
parameters used in this figure are the same as in Fig. ||. Here, we also defined that is real, and 
hence the valleys coming from the Hubble A-term lie along (arg(</>) = vr/4, 37r/4, 5tt/4, 7tt/A). We 
see that the phase of the (j> field starts to oscillate around the bottom of the valley at H ~ Hph- 
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FIG. 4. The evolution of the baryon asymmetry — — = — estimated by a numerical cal- 

s 23 s 

culation. Parameters used in this figure are the same as in Fig. and pi Here, we also take 
Tr = 10 4 GeV. The thick line denotes the estimated value using Eq. (|3l|). We can see that the 
asymmetry begins to oscillate with constant amplitude at H ~ H p h and it is fixed when the <p field 
starts the oscillations around the origin. The analytical estimation of the baryon asymmetry agrees 
well with the obtained value by the numerical calculation. 
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FIG. 5. The plot of the parameter regions which are free from the early oscillation in m u i — Tr 
plane. The regions below the four solid lines are free from the early oscillation. These four 
lines correspond to the breaking scale of the U{1)b-l-> v = 10 16 , 10 15 , 10 14 , 10 13 (GeV) from 
left to right, respectively. The shaded region corresponds to the present baryon asymmetry, 
ns/s ~ (0.4 — 1) x 10~ 10 . We have taken S e ff, |c^| = 1 in this figure. Here, we assume the 
existence of or the /x-term when the amplitude of the <f> field is of the order of the U(1)b-l 
breaking scale v. From this figure, we see that if we take the breaking scale of the U(1)b-l, 
v ^ 10 14 GeV, the early oscillation of the cf> field H osc > H p h can be avoided in most of the regions 
of parameter space. 
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